Hysteresis in the electronic transport of V2O3 thin films: non-exponential kinetics and 

range scale of phase coexistence. 
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The thermal hysteresis of the electronic transport properties were studied for V2O3 thin films. 
The temporal evolution of the resistance shows the out-of-equilibrium nature of this hysteresis with a 
very slow relaxation. Partial cycles reveal not only a behavior consistent with phase coexistence, but 
also the presence of spinodal temperatures which are largely separated. The temperature spreading 
of phase coexistence is consistent with the bulk phase diagram in the pressure-temperature plane, 
confirming that the film is effectively under an effective pressure induced by the substrate. 
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INTRODUCTION 



The physical properties of V2O3, cited as a prototype 
of strongly correlated systems, have been extensively 
studied especially for its spectacular metal-to-insulator 
(M-I) transition [l|, @] • Experimental evidences of a 
strong pressure dependence of physical properties have 
been also accumulated. Consequently, the phase diagram 
in the pressure-temperature (P,T) plane was established 
[3| . A lot of work has been devoted to the M-I Mott tran- 
sition which can be observed when V2O3 is subjected to 
a chemical pressure by Chromium substitution [5j . Such 
a negative pressure can be approached in thin films when 
growing V2O3 under tensile stress on LiTa03 substrate 
[6|. The other M-I transition is associated to a change 
from a rhombohcdral symmetry to a monoclinic symme- 
try and occurs at T ~ 160 K under atmospheric pressure 
[7( . This structural transition leads to an increase of the 
unit-cell volume with range going from 1% to about 3.5% 
Q. This transition has a strong first order character, and 
a 10 K thermal hysteresis is often reported. When a hy- 
drostatic (or uniaxial) pressure is applied, the transition 
temperature decreases, and above a critical pressure of 
roughly 26 KBars the metallic state is completely stabi- 
lized down to the lowest temperature 0]. As a conse- 
quence, this strong pressure-dependence of the M-I tran- 
sition is a motivation for thin film studies since large 
effective pressure can be expected by substrate-induced 
strains [9( . Recently, thin films of V2 O3 were grown on 
sapphire substrates, using pulsed laser deposition tech- 
nique Under certain growth conditions, the M-I 
transition of the films can be suppressed at the macro- 
scopic scale for peculiar thicknesses [ijj]- Instead of the 
insulating- like resistance for T < 160 K, a broad ther- 
mal hysteresis surrounding a non-monotonic (dR/dT) is 
indeed observed. To explain the (dR/dT) shape, we 
demonstrated that below 160 K, the samples exhibit an 
electronic phase separation between metallic and insulat- 
ing mesoscopic zones whose concentration ratio depends 
on the thickness [3| • Indications of a broad distribution 
of critical temperatures have also been highlighted and 



could be a cause for a large hysteresis. It is thus nec- 
essary to understand this distribution, and consequently 
its role in the shape of the hysteresis. In this paper, we 
show using the technique of minor loops that this broad 
phase coexistence is responsible for the hysteresis and for 
non equilibrium transport properties over a large temper- 
ature scale. 



EXPERIMENTAL 

Briefly, a series of V2O3 thin films, with various thick- 
nesses (i), was epitaxially deposited on (OOOl)-sapphire 
substrate using the pulsed laser deposition technique. 
Details of the growth can be found elsewhere [IH . Usu- 
ally, thin films are under an effective pressure below a 
critical thickness, where the substrate causes the film to 
grow with uniform elastic strain [ijj]. In the case of a 
large mismatch between film and substrate in-plane lat- 
tice parameters, the strain is quickly released in the film, 
leading to the so-called misfit dislocation regime. No 
critical thickness is then observed down to the ultra-thin 
film limit. It is the case here where the lattice mismatch 
between V2O3 and sapphire is about 4 %, as confirmed 
by X-ray diffraction. The tendency is even to recover the 
bulk parameters at low thicknesses [3]. 

The electronic resistivity (p) is measured using the 
four probe geometry in a commercial PPMS (Physical 
Properties Measurement System). The resistivity was 
estimated with the thin rectangular slice approximation: 
p = Rir/ln2.tRi(W/L), where R is the resistance, W the 
width, L the length and Ri(W/L) is the correction due 
to the finite siz e 11511 . A microbridge, whose preparation 
is described in [1 41 ] . from a llnm thick film with a size 
of 100 x 300pm 2 was also measured. The detailed anal- 
ysis of partial cycles and of resistance relaxation were 
performed in this microbridge. For large current, the 
voltage-current characteristics show non linear behavior, 
as sometimes observed in samples with phase separation 
[l6j |. However, we limit our study well inside the ohmic 
regime which is within the range / = 1-500 pA. We also 
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Figure 1: Temperature-dependence of the normalized re- 
sistance measured during a thermal cycle (cooling-heating, 
rate=0.5 K/min) for 3 typical thicknesses of V2O3 films. A 
large thermal hysteresis is observed from 20 K < T < 150 K. 



Figure 2: Evolution of the thermal hysteresis width as a func- 
tion of the V2O3 film thickness. The constant trend confirms 
that no size effect is observed. 



checked that the relaxation of the resistance was uncorre- 
cted with the temperature-dependence of the resistivity 
and that Joule heating can be avoided. The results pre- 
sented hereafter were measured for a typical current of 
IfiA. 



FIRST-ORDER TRANSITION IN V 2 3 THIN 
FILMS 

Recently, we have shown a strong thickness- 
dependence in the properties of the V2O3 thin films 
[121] . In particular, we have observed that the M-I transi- 
tion can be tuned by changing the thickness of the films 
from 2 nm (±1 nm) up to 300 nm (±5nm). For exam- 
ple, the M-I transition is observed for ultrathin samples 
(t < 10 nm) with lattice parameters close to the bulk 
values, whereas thicker samples (t > 22 nm) present a 
quasi-metallic behavior down to low temperatures. For 
intermediate thicknesses, the situation is more complex 
: the films exhibit an intermediate behavior with a clear 
maximum in the R{T). All samples present an thermal 
hysteresis AT in the range 20K-150K, but the magni- 
tude of the hysteresis is maximum for the intermediate 
thicknesses. These various behaviors are illustrated in 
figure [T] where R(T) curves are plotted for three typical 
thicknesses. 

The M-I transition in V2O3 single crystals occurs with 
strong first order characteristics Generally, the 

hysteresis of a first order phase transition results from 
the possibility of supercooling (superheating) metastable 
states [18(- The hysteresis observed in our films could 



indicate such a first order transition, but its width is 
surprisingly large (AT « 130 K compared to 10 K for the 
bulk samples) . The broadening of the hysteresis seems to 
be specific to the thin film form, and may originate from 
several factors as described below. The first factor is that 
the decreasing of the system dimensions (or system thick- 
ness) can induce a finite size effect, which can lead to the 
so-called "rounding field" As a consequence, the 

width of the first order transition is strongly dependent 
on the thickness, and AT C cx t -1 . This has been indeed 
used to explain the thermal broadening of the M-I tran- 
sition in VO2 films 20]. In our V2O3 samples, we do 
not observe any significant thickness-dependence of the 
thermal hysteresis width as seen in figure [2] The second 
possibility is that the growth generates quenched disor- 
der, which can modify a first order transition to a con- 
tinuous one. Actually, this has been shown for a dimen- 
sionality D = 2 on the basis of general grounds [2l| . For 
D = 3, the quenched disorder can also smooth the tran- 
sition even if this latter still retains first order character- 
istics. To understand the role of such disorder upon the 
width AT of the hysteresis, our samples were annealed 
under argon atmosphere, since such thermal treatment is 
usually used to order a sample (as in manganites [13]). It 
is unlikely that the quenched disorder (cationic or oxy- 
gen disorder) plays a strong role in our films because, 
even if the residual resistivity slightly changed by about 
10 %, the shape of the R{T) and the width AT of the 
hysteresis remain unchanged after thermal treatement. 
Finally, a third source of heterogeneity in our films can 
be the clamping of the unit cell that can induce (likely 
heterogeneous) strains [23]. Using a two-parallel resis- 
tor model, we explained that the peculiar shapes of the 



3 




t (sec) t (sec) 



Figure 3: Normalized resistance measured as a function of 
time at different temperatures obtained during a cooling on 
the microbridge. The time dependence evidences a metastable 
character. 

R(T) curves (or the non-monotonic thermal variation of 
the resistance slope dR/dT) obtained when cooling the 
sample arises from the coexistence of metallic and insu- 
lating zones in the sample However, the link be- 
tween this phase coexistence and the large hysteresis re- 
mains unexplored. The large temperature-range of phase 
coexistence can indicate a distribution of critical tem- 
peratures at a macroscopic scale. Since the reduction of 
critical temperature is directly correlated with the pres- 
sure and considering that the substrate acts as a strain 
field, this means that some areas of the film are strongly 
strained, whereas other ones are partially relaxed. In 
other words, there should be locally a distribution of the 
pressure within the films. We will now focus on the na- 
ture of the thermal resistance hysteresis and show how it 
can be correlated on the (P,T) phase diagram. 

TIME-DEPENDENT EFFECTS AND NON 
EXPONENTIAL KINETICS 

To understand the nature of the hysteresis, we have 
measured on a microbridge of 11 nm thick film the time- 
dependence of resistance. A slow relaxation can be ob- 
served for 20 K < T < 150 K (Fig. [3]). The sign of 
the relaxation depends on whether the metastable state 
is stabilized after heating or cooling the sample. 

To fit the time-dependence resistance (R(t)), different 
functional forms can be proposed : exponential, stretched 
exponential or logarithmic time dependence. Here, the 
relaxation is clearly non-exponential. The simplest expla- 
nation of such non-exponential kinetics, which is ubiqui- 
tous in disordered systems, is that the relaxation can be 
intrinsically exponential in different domains, but that 



Figure 4: Evidence of the non-exponential kinetics of the re- 
sistance relaxation at T=80K (The total acquisition time of 
about one day). The data are presented in a semi log-scale 
to emphasize the logarithm relaxation. The solid line is a fit 
R(t)/R(t = 0) = 1 + A.ln(t + T ) (A = 0.0004,r = 250). In 
the inset are shown the stretched exponential fit (dashed line) 
and the logarithm fit of the low time data. 



the relaxation time is not unique at the sample scale. 
This is in agreement with the idea that different do- 
mains with different critical temperatures are present in 
the sample. For non interacting domains, a large statisti- 
cal averaging of exponential relaxation times leads to the 
stretched exponential in the form R(t) oc exp(±t/ro)P 
(where tq is a relaxation time and j3 < 1 the relaxation 
exponent) Q • If the acquisition time is moderate (typ- 
ically of the order of one hour), then such a stretched 
exponential provides a good fitting of the date. The non- 
exponential time dependence obtained ((3 « 0.2) confirms 
a non Debye relaxation, and indicates that the distribu- 
tion function of characteristics time is rather flat. The 
second case to describe the data leads to a logarithm 
relaxation of the form R(t) oc ln(t + tq). In order to 
distinguish between the two fits, very long relaxation 
times (t ~ 1 Day) were performed, and the resulting 
curves are presented in figure [4] The plots show that 
the relaxation is even more slower than expected from 
the stretched exponential. The plot clearly evidences a 
logarithmic-dependence of the resistivity as a function 
of time, whereas the stretched exponential is not sat- 
isfactory for long periods (unless we do not add a sec- 
ond stretched exponential somewhat arbitrarily to fit the 
data) . 

Logarithm relaxation is usually observed in complex 
and strongly interacting materials [24| . For example, it 
is observed in glasses, whether they are structural, mag- 
netic or electronic, and this relaxation corresponds to 
the ageing of the system. To check the hypothesis of 
glassy properties in our sample, we tried to observe spe- 
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cific properties of glassy states such as memory and re- 
juvenation [25| . Note that in the case of a glassy state, a 
memory effect is usually seen after a negative tempera- 
ture cycling, because the spent time at the low tempera- 
ture does not influence on the ageing of the system. Here, 
a negative temperature cycling during the relaxation (80 
K — ► 75 K — ► 80 K) stops the resistance relaxation and 
then helps the system to reach its equilibrium. Since we 
do not observe genuine glassy characteristics, the slow 
kinetics are not due to the fact that the system always 
explores metastable states. The relaxation is rather due 
to the nucleation of M-I domains with a large number 
of transition temperatures into a M matrix (when cool- 
ing the sample). This broad distribution of transition 
temperatures is the reason for the large number of relax- 
ation times inferred by the logarithmic relaxation. This 
relaxation is very slow, showing that the kinetics are ex- 
tremely constraint. Wc think that this is due to the long 
range strain effect induced by the substrate because of 
the clamping of the unit cell at the bulk transition tem- 
perature. This kind of relaxation is similar to the kinetics 
of phase transformations subjected to pressure [2(| . 



MINOR LOOPS AND PHASE COEXISTENCE 

It is now necessary to determine the region of phase co- 
existence and the stability limit of the metastable states, 
e.g. the spinodal limit. A way to confirm and to get 
more informations on the phase coexistence is to per- 
form subloops in the hysteresis cycle. The reason is the 
following. As far as the transformation from one state to 
another is not complete, partial hysterctic cycles can be 
observed. 

Typical partial cycles measured on the microbridge are 
shown in figure [5] and figure [Bj. For the cycling sequence 
high T— dow T— >high T shown in figure El we observe 
that the metastability begins at T as 150 K and ends 
at T as 20 K. However, the reverse cycling sequence low 
T— »high T^low T shown in figure [5] exhibits a quite dif- 
ferent behavior. The three first curves are fully reversible 
(label 1, 2 and 3 in the figure [6]) whereas the 70A"-one 
(label 4) shows irreversibility. Since it is necessary to 
heat at temperature T > 70 K to observe that a ther- 
mal hysteresis begins to develop, this means that phase 
coexistence begins at T as 70 K and ends at T as 150 K. 
This asymmetry appearing for partial cycles and depend- 
ing on the initial temperature indicates that below 70 K, 
the high resistance state is strongly favored and can be 
regarded as the stable state (even if some parts of the 
sample can be stabilized in the low resistance state by 
supercooling from high temperature). Thus, the analysis 
shows that the zone of phase coexistence extends only 
from 70 K to 150 K. From the (P,T) bulk phase diagram 
of the V2O3, the transition line which determines the M-I 
transition can be described as a straight line with a slope 
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Figure 5: Partial cycles of resistance versus temperature with 
the sequence high T-low T-high T measured on the micro- 
bridge. See text for details. 
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Figure 6: Left: Partial cycles of resistance versus temperature 
with the sequence low T-high T-low T on the microbridge. 
Right: Zoom of the low temperature partial cycles. The re- 
turn in the cycles is made at the temperature ?X)K (label 1), 
45A (label 2), 50 A" (label 3), and 70K (label 4) (see text for 
details). 



AT C /AP as -3.5K/Kbars for T > 70K and a quasi- 
horizontal line with a slope AT C /AP as ~AQK / Kbars for 
T < 70 K 8] (FigfTJ). The pressure-induced at the tran- 
sition has been thus calculated to be AP max as 27KBars 
for a complete blocking of the transition [l4[ . This value 
is higher than the critical pressure (26KBars) reported 
in V2O3 single crystals by McWhan et al § and shows 
that a fully strained film should not exhibit any M-I tran- 
sition. Since a first order transition leads to an hetero- 
geneous state when the pressure increases at a constant 
volume, we assume that our film sustains a very broad 
and continuous effective pressure distribution (from P as 
to P as 27KBars, where each pressure is statistically 
equivalent). Therefore, the film should also display a very 
broad critical temperature distribution ranging from K 
to 150 K. However, from the (P,T) diagram, a simple 
estimation shows that the phases whose critical temper- 
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Figure 7: Relationship between pressure-temperature (P,T) 
plane of V2O3 bulk material extracted from [i| and a his- 
togram of the domain proportion exhibiting one critical tem- 
perature (see the text for details). 



atures are in the range < T c < 70 K are statistically 
unfavorable. As a consequence, the temperature range 
where the phase coexistence occurs is going from 70 K 
to 150 K, in agreement with the analysis of the partial 
cycles. This is summarized in figure [7j Even if partial 
cycles were not performed for all the thicknesses investi- 
gated, the hysteresis cycles were always similar in terms 
of characteristic temperatures (see Fig (2]). This suggests 
that our analysis is not thickness-dependent as far as the 
width of the phase coexistence is concerned. Finally, one 
has to note that V2O3 crystals form cracks because of 
fractures, which occur during the volume change at the 
M-I transition. This is not the case in our films. In bulk 
V2O3, M and I phases correspond to different structural 
symmetries, which transform into each other via a (local) 
first order transition. Some indications of the thermoe- 
lastic martensitic effect importance in the M-I transition 
in V2O3 have been reported using acoustic emission mea- 
surements [27]]. We propose that in the particular case 
of V2O3 films, substrate enables a releasing of the strain 
by inducing a heterogeneous state, as it is the case in a 
martensitic scenario, where a large number of small do- 
mains increases the spread in transition temperature of 
each domain 12811. 



CONCLUSION 

We have presented a detailed study of thermal hys- 
teresis present in the transport properties of V2O3 thin 
films. Non-exponential relaxation and partial cycles are 
observed, showing a broad region of phase coexistence. 
Using partial hysteresis cycles, this zone of the phase 
coexistence has been precisely established. The limit of 
metastability corresponds well to what can be expected if 



the film is under an effective pressure induced by the sub- 
strate. We believe that such results will be of the general 
interest for the thin film physics that display metal-to- 
insulator transition. 

This work was carried out in the frame of the Eu- 
ropean project CoMePhS (NMP4-CT-2005-517039) sup- 
ported by the European Community and by the CNRS, 
France. 
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